Charge transfer statistics of a molecular quantum dot 
with strong electron-phonon interaction 
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We analyze the nonequilibrium transport properties of a quantum dot with a harmonic degree of 
freedom (Holstein phonon) coupled to metallic leads, and derive its full counting statistics (FCS). 
Using the Lang-Firsov (polaron) transformation, we construct a diagrammatic scheme to calculate 
the cumulant generating function. The electron-phonon interaction is taken into account exactly, 
and the employed approximation represents a summation of a diagram subset with respect to the 
tunneling amplitude. By comparison to Monte Carlo data the formalism is shown to capture the 
basic properties of the strong coupling regime. 
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In the past decades, the miniaturization of electric cir- 
cuits has crossed the divide between the microscale and 
the nanoscale. Molecular and atomic electronics are no 
longer mere theoretical concepts. A wide variety of ex- 
perimental setups has been developed for the exploration 
of the electronic and mechanical properties of nanometer- 
sized objects like carbon nanotubes, Cgo-fullerenes and 
complex molecules. It has become possible to connect 
these samples to mesoscopic environments and to inves- 
tigate their transport properties^— 

The electronic structure of these nanosized objects is 
best captured by the concept of a quantum dot, i.e., by 
an arrangement of energy levels which correspond to the 
molecular orbitals of the device. One of the most promi- 
nent and fundamental models for the theoretical descrip- 
tion of quantum dots is the Anderson impurity model, 
which accounts for the tunnel coupling to noninteracting 
electron reservoirs and for the local Coulomb interaction 
between the electrons populating the quantum dot*£ In 
the case of contacted molecules, where charging is often 
accompanied by structural deformations of the molecule 
itself, however, this model is often an oversimplification. 
For a more realistic description, an explicit consideration 
of the coupling to vibrational degrees of freedom is nec- 
essary. This is accomplished by the Anderson-Holstein 
model (AHM)J^ 

In its full extent, the AHM captures a huge variety of 
physical phenomena. Its physical properties depend on 
several energy scales, e.g., temperature, charging energy, 
hybridization energy, level spacing and electron-phonon 
interaction strength. These define many interesting and 
physically distinct regimes in parameter space. In this 
paper we are mainly interested in the effect of electron- 
phonon interactions on the charge transport through a 
contacted molecule. The model can therefore be simpli- 
fied to contain a single electronic level (thus neglecting 
the spin degree of freedom as well as the charging en- 
ergy) linearly coupled to a local (Holstein) phonon, i.e., 
a bosonic oscillator degree of freedom with a single fre- 
quency. Even this simplified model, which in the fol- 



lowing will be referred to as AHM, offers rich physics. 
Whereas the conductance and the nonlinear I — V- 
characteristic of such a system can be approached by a 
number of methods, such as diagrammatic Monte Carlo 
schemes,— rate equations , 11 : 16 perturbation theory,— ~— 
and P(E) theory— its full counting statistics (FCS) is 
well understood only in the limit of weak electron-phonon 
coupling 

In this paper we would like to extend these results 
and present a calculation of the FCS beyond the weak- 
coupling limit. One possible experimental setup in 
which strong electron-phonon interaction can be reached 
is a quantum dot embedded in a suspended carbon 
nanotube ) 11 ' 29 : 30 as depicted in Fig. [T] The electrodes 
as well as the quantum dot are made of a single car- 
bon nanotube subject to a bias voltage V. The elec- 
tronic level structure of the quantum dot can be tuned 
by an additional backgate. Because of its simple struc- 
ture, the vibrational modes of such a quantum dot are 
well understood^ However, the model employed is fairly 
general. Depending on the parameter regime, it also al- 
lows the desription of transport through molecules con- 
tacted using mechanically-controlled break junctions 2 * 9 - 
and STM tips^* 3 - 2 - as well as in nanoelectromechanical 
setups^ 3 - 

The starting point of our calculations is the Hamilto- 
nian of the AHM, 



H — Hq + Ht + H e i. p b- 



(1) 



Hq is the Hamiltonian for the uncoupled degrees of free- 
dom of the electrodes, the single-level quantum dot, and 
the bosonic mode, 



H = H L [ip L ] + H R [4, R ] + e dU + ntfb. 



(2) 



The first two terms describes the electrodes in the lan- 
guage of noninteracting electron field operators ipL,R. (x) 
which are held at chemical potentials hl,r- This is 
achieved by a bias voltage V = hl^^r applied across the 
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FIG. 1: (Color online) Sketch of a suspended carbon nanotubc 
quantum dot. The fields i/jl,r{x) describe semi-infinite leads, 
and d is the electron annihilation operator on the quantum 
dot. The gate voltage Vg allows for an adjustment of the 
energy levels of the dot. The bias voltage V induces a finite 
charge current across the dot. P (Q) denotes a measurement 
apparatus for the charge transfer statistics. 



contacts (we use natural units e = h = fee = 1 through- 
out). In the case of noninteracting leads, the only prop- 
erty of Hl,r which is necessary for further calculations 
is the local tunneling density of states p(oj). As a simpli- 
fication, we work in the wide fiat band limit p(to) = po, 
which corresponds to free fermions iI>l,r with linear dis- 
persion and infinite bandwidth, but any other shape of 
p{uj) can be treated in the same way. The third term 
in Ho comes from the single electronic level of the quan- 
tum dot with bare energy eo and associated creation and 
annihilation operators and d. The last term in Hq 
describes a phonon mode with frequency fl and bosonic 
creation and annihilation operators and b. 

The tunneling Hamiltonian Ht describes the hopping 
of electrons from the leads to the dot and back: 



H ^ = i E [v4(* = o)d+h.c], (3) 

m—L.R 



where 7 is the tunneling amplitude and we assumed 
symmetric coupling to both leads. An asymmetric cou- 
pling requires only small adjustments in the calculations 
but significantly complicates notation. The last term in 
Eq. ([I]) is the electron-phonon interaction, 



H c _ ph = gdU (fct + b) 



It linearly couples the dot occupation operator d^d to the 
displacement of the harmonic oscillator Q ~ fot _|_ 5 with 
a coupling strength g. 

The fundamental quantity describing low-frequency 
electronic transport is the probability distribution func- 
tion P (Q) of transferring Q units of charge during a mea- 
surement time T, which we assume to be the longest time 
scale. Physical observables can then be calculated as av- 
erages with respect to this distribution function. The 
expectation value (Q) is related to the average current, 
(Q) = (/) T. The second cumulant ((Q 2 )) = (Q 2 ) - (Q) 2 
is directly related to the noise power ((Q 2 )) = ST, where 
S = ±J du((I(t)I(0) + I(0)I(t))). Instead of calculat- 
ing the probability distribution function itself it is of- 
ten more convenient to calculate its cumulant generat- 
ing function (CGF) In* (A) = hi^g c ix ®P (Q). It has 



been shown^^ that In x(A) can be expressed in terms of 
Kcldysh Green's functions (GFs) and this formalism has 
been successfully applied to a wide range of transport 
problems. The fundamental expression for calculating 
the CGF isM 



x (A) = (T c e- i /cd*T A (*)^ 



(5) 



where Tq denotes time-ordering along the Keldysh con- 
tour C and the expectation value is written in the interac- 
tion picture with respect to the Hamiltonian Hq + -ff c i-ph- 
The tunneling operator T\ is given by 



T - 7 
J* — — 



iA/ Vl(0R 



■ u/ Vi(0)d- 



h.c. 



(0) 



The counting field A is explicitly time-dependent on the 
Keldysh contour: A(t € C±) = A . At the end of the 
calculation one has to replace A + = 
shown in Ref. (37| . in the limit T 
expression holds 



A = A. As was 
00 the following 
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The A-dependent expectation value is defined as 
'T c ■■■ e-i/cd* 1 ^ 



X(A+,A") 



(7) 



(8) 



We proceed by defining the exact A-dependent dot GF 
D\ and the local (taken at x = 0) free electrode GF g m 
(where m = L, R) in Keldysh space by 



D x (t,t') = -i(T c d(t)d\t')) x7 

g m (t, t 1 ) = -i (T c ^ m (0, t)i>l(0, t')) Q , 



(9) 



Then, Eq. ([7]) can be expressed as the following convolu- 
tion, 



91nx(A",A+) iT7 2 



(lid 

2^ 



(10) 



E 



c-^-+h.9+-h 



■(w)£>t-(w) 



where A = (A - — A + )/4. The usual route for the 
limit of weak electron-phonon coupling would be an 
expansion of the dot GF in the interaction strength 
g. Alternatively, the strong coupling regime can con- 
veniently be approached using a Lang-Firsov (polaron) 
transformation^, JJ = exp \ad)d(b^ — b)\ . One can 
check that the Hamiltonian transforms to U(Hq +T\ + 
H eUph )tf =H' a +T' x + H' cl _ ph , where 



H'o 



Ti = 



7 



H R + (e - 
-aftdU(b 



E 



,-imA/ 



m=L,R.= -,+ 

gdU (b + tf) - 



v /VL(0)e Q(6t ~ b) d + h, 
2agd)d. 



11) 
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FIG. 2: Feynman diagrams for the contributions = -Do, 
D\ and D\. The solid lines represent the free dot GF d(t,t') 
and the dashed lines the free lead GF gL,R.(t,t'). The internal 
time variables U are integrated over while t, t' (denoted by the 
filled circles) are external time variables. 

Setting a = g/£l removes the electron-phonon interaction 
and leads to the polaron shift of the dot energy, eo — > 
eo — g 2 /n. The polaron shift can be taken care of by 
appropriate gating and, for ease of notation, we use eo 
to denote the shifted dot level from now on. Moreover, 
the transformed tunneling operator now contains the 
dressed dot operator 

D:=e 9 ^- bS >' n d=:Xd. (12) 

The operator creates a particle on the dot which is 
dressed by a phonon cloud X . The next step is to ex- 
press the dot GF in the language of the new transformed 
Hamiltonianj 19 ' 39 One finds 

D x (t,t') = -i(T c D(t)D\t')) x 

= -\(T c d{t)d\t')X{t)X- 1 {t')) x , (13) 

where the A-expectation value now has to be calculated 
using the rotated Hamiltonians H' and T' x . In the next 
step we perform a formal expansion in the tunneling 
Hamiltonian T x and find 

oo 

D x (t,t') = J2l 2n DT(t,t') , (14) 

n=0 



where D^{t, t') is the sum of all terms to 2roth order in 
the tunneling amplitude. As an example, the diagram- 
matic structure of the contributions D x , D\ and D\ is 
shown in Fig. [5] In addition to the free dot GF, 

d(t,t') = -i(T c d(t)J(ti)) , (15) 

the calculation of the individual orders in- 
volves correlators of the form A2„(<i, . . . , t^n) = 
(T c A(t 1 )X- 1 (t 2 )...A(i 2n _ 1 )X- 1 (i 2 , l )>o. Since the 
A-operators are exponentials of free boson operators, the 
calculation of this correlation function is straightforward: 

2n 

A 2 n(ti,...,t 2n ) = Y[A(t t -t J ) , (16) 

i<j 

where A(t — t') is a function in Keldysh space, 
A(t-f) = {A«(t-f)} fcil=± 

f K (\t-t'\) K{e -t) \ > ( 17 ) 

-\*{t-t) K(-\t-t'\)J 

and n(t) is defined as 

«(t) = cxp {-a 2 [(e int - 1) n B + (e" int - l) (n B + !)]}• 

(18) 

The uncoupled phonon occupation number ne = (b^b)o 
accounts for the initial occupation of the harmonic os- 
cillator states. If the oscillator is coupled to a thermal 
environment^ e.g., to a substrate or a backgate, tib is 
a temperature-dependent distribution function. In this 
case, it is often sensible to assume that it given by an 
equilibrium Bose distribution, iib(T) = (c n ^ T — l) -1 , 
where T is the bath temperature. In the limit T — > iib 
approaches zero. Note that tib denotes the phonon num- 
ber in the absence of coupling to the dot. It is generally 
different from the highly nontrivial, nonequilibrium occu- 
pation number that emerges as a consequence of the cou- 
pling to the dot, and which is governed by the transport 
processes in the lead-coupled systemj 28 i 40 The function 
n(t) can be expressed as a Fourier series 



e _ Q 2 [2 „ B(T) + 1] ^ (nB + 1) 



n— — oo 



T = 0. 

e nO/2T e -inf2t J 1 > 



(19) 



where I n denotes the nth order modified Bessel function. 
There is a crucial difference between the T = and T > 
expansion. In the former case, only positive phonon 
numbers n occur. This is natural because in this case, the 



phonons can only be excited. At finite T, in contrast, n 
runs over positive and negative integers because now the 
phonons can be emitted or absorbed. In the limit T — > 
both expressions have to coincide. This can be verified 
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FIG. 3: Linear diagram resummation scheme. The single solid 
lines denote the free dot GF d(t,t') and the dashed lines de- 
note the free electrode GFs gL,n(t,t'). Wiggly lines represent 
the phonon cloud propagator A(t,t'). This approximation al- 
lows to derive the Dyson equation in the second line, where 
Ql.r acts as self-energy. 



by considering the x —> limit of the modified Besscl 
function I n (x) = x n { [2'T(n + 1)]" 1 +0(x 2 )} where T(n) 
denotes the gamma function^ 

So far, no approximations have been performed. How- 
ever, the complicated structure of the diagrams (see 
Fig. [2]) makes an exact solution of the problem diffi- 
cult. Moreover, a simple perturbative expansion in the 
tunneling amplitude 7 is inconvenient as the individual 
contributions are highly divergent. One of possible reg- 
ularization methods is to couple the quantum dot to a 
fictitious bath to enforce a hybridization of the sharp dot 
level (aka Wigner-Weisskopf regularization)4i We take 
another path and calculate the dot GF using a linear ap- 
proximation of the diagrammatic series (see Fig. [3]). It 
allows us to derive the following Dyson equation in fre- 
quency domain, 



D 
D 



■7 D 
D SD^ 



e 2 gR 



Dr 



(20) 



Here, we used the conventional matrix notation in 
Keldysh space, [Do(w)]h = Dq 1 {lo) for k,l = ±, and 
similar for the self-energy S. The function Dq 1 (lu) is 
the convolution of A kl (uj) with the free dot propagator 
d kt (oj), i.e., D% l {t,t') = A kl (t,t')d kl (t,t'). In this ap- 
proximation the dot GF can be determined exactly. A 
lengthy but straightforward calculation reveals that the 
CGF has the form of the Levitov-Lesovik formula^ 



lnx(A) =rj ^lnjl + TH 

x [n L (1 - n R ) e iA + n R (1 - n L ) e" iA ] 




FIG. 4: (Color online) Transmission coefficient T(lo). Main 
graph: Transmission coefficient for T — eo = 0, fi/T = 5, 
and g/T = 1, 5, 10 (solid, dashed, dashed-dotted line, respec- 
tively). Inset: Transmission coefficient for fixed g/T = 1, 
Q/r = 5 and eo = for varying temperatures T/T = 1, 10 
solid, dashed line. 



pansions in the regimes T — and T > 0. 

2 . . 2n -1 

„ — a a 1 



f 



%\ oj — £{)—nQ, 
-a 2 (2„ B + l) g /42aV« B (n B + l)1e 



Tif2/2T 



n— — 00 



oj — eQ—nQ 



T = 0, 

T > 0. 

(23) 



Again, in the case T = 0, only processes involving the 
emission of phonons are allowed. In contrast, for T > 0, 
thermally excited phonons can be absorbed. The trans- 
mission coefficient T(lu) is made up of a sequence of peaks 
(see Fig. gj at the energies Nfi for T = or ZVL for 
T > and is properly normalized, J duJ \u) / V = 1. By 
assuming a Lorentzian shape, the width of the peak at 
lu = nVt can roughly be estimated as 2Te~ a a 2n jn\ and 
2re -a 2 (2r lB +i) /n [ 2a 2^ nB (n B + 1)], for T = and T > 
respectively. 

Interestingly, both perfect and zero transmission are 
possible in the system. This is due to the special struc- 
ture of T((jj), which is equivalent to the transmission 
coefficient of a system with an infinite number of spin- 
degenerate dots at energies eo + nil coupled to t he lead s 



(21) 



HH). 



in parallel (for the double dot system see, e.g., 
This occurs as a consequence of the linear approximation: 
every single electron tunneling through the system takes 
along its polaron cloud, leaving the dot in exactly the 
same state as before the tunneling event. That means 
that the resonance condition is given by u> = eo + nO. 
The antiresonance (complete transmission suppression) 
emerges as an interference effect in precisely the same 
way as in the double-dot setupi 42 i 43 The above physical 
picture implies that the dwell time of electrons on the dot, 
which is on the order of T -1 , must be long compared to 
the inverse of the phonon (de)excitation rate. Since the 
with r = 27rpo7 2 - The function /(w) has different ex- latter is roughly proportional to g, we expect g 3> T to be 



where nL,n(u) are Fermi distribution functions in the 
lcft(right) lead with chemical potentials /j,l,r = ±V/2 
and T(u) is the effective transmission coefficient 



r 2 



r 2 



(22) 
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FIG. 5: (Color online) Current in single particle approxima- 
tion, linear approximation and Monte-Carlo (MC) simulation 
data of Ref. [lj] for T/T = 0.2. The dashed lines represent 
the single particle approximation as proposed in Ref. [10] , the 
solid line is the linear approximation and the dots represent 
the diagMC data. The energy level of the dot is held at eo = 0. 
The assignment of the data (from top to bottom) is the fol- 
lowing (g/T = 2, fi/r = 5, red), (g/T = 4, ft/T = 5, green) 
and (g/T = 4, fi/r — 3, blue). The systematic deviation 
of the analytical results from the diagMC data for large V/T 
is due to finite bandwidth necessary to carry out numerical 
simulations. 

a necessary requirement for the validity of our approxi- 
mation. This indeed implies a strong electron-phonon 
coupling. 

At all temperatures one observes an exponential sup- 
pression of the peak width with the coupling a 2 = 
(g/fl) 2 . This leads to the well known Franck-Condon 
blockade, where sequential tunneling is exponentially 
suppressed and tunneling accompanied by phonon ab- 
sorption/emission is preferred^ This is also observable 
in the current. In Fig. [5l the current is depicted for 
different electron-phonon coupling constants g. For in- 
creasing coupling the step heights modifies nonuniformly 
and transport through states with higher n— phonon ex- 
citation, n « (g/fi) 2 , gets more pronounced. The ef- 
fect of temperature is similar to the phonon coupling 
strength a (this is obvious, because a 2 (2hb + 1) or 
2a 2 ?iB (pb + 1) always acts as an effective coupling 
strength). In order to assess the quality of our approxi- 
mation we compared the calculated I — V characteristic 
with the one from diagrammatic MC data of Ref. flE] ]. 
In the regime of small to moderate V and O, as well as 
for g > r, our scheme indeed turns out to yield a better 
approximation than the single-particle approximation.— 

The noise is plotted in Fig. [5] Similar to the I—V char- 



acteristic, is also shows a step-like behavior. For strong 
coupling there is an additional feature: in the steps we 
observe an additional plateau. For large voltages the 
noise approaches the usual unitary limit, in our units 
S = I. Unsurprisingly, our approximation, being valid 
for not too large currents, shows no sign of enhanced noise 
due to the predicted avalanche-like transport behavior.— 
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FIG. 6: (Color online) Main graph: the noise in the linear 
approximation is depicted at T = for g/T = 1, 2, 4, 6 (solid, 
dashed, dashed-dotted, dotted line, respectively). The re- 
maining parameters are set to Q/T — 5 and eo = 0. In- 
set: the noise for temperatures T/T = 1, 5, 25 (solid, dashed, 
dashed-dotted , respectively). The remaining parameters are 
g/T = 1, fi/r = 5 and e = 0. 



To include higher-order correlations (see for example the 
nearest-neighbor crossing approximation for the current 
in Ref. [HI) would be a task for the future. 

In conclusion, we developed an approach to calculate 
the FCS of the Holstein polaron dot in a strong coupling 
regime. Using a linear approximation, we derived an ana- 
lytical Levitov-Lcsovik formula for the cumulant generat- 
ing function with an effective, properly normalized trans- 
mission coefficient. Our approach yields predictions for 
zero temperature as well as for finite temperature, where 
the phonon is assumed to be thermally equilibrated. 
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